Partial magnetization plateau emerging from a quantum spin ice state in Tb2Ti207 
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The quantum spin ice model applied to Tb2Ti2 07 predicts that magnetic fields applied along the 
[111] axis will induce a partial magnetization plateau [H. R. Molavian and M. J. P. Gingras, J. Phys.: 
Condens. Matter 21, 172201 (2009)]. We test this hypothesis using ac magnetic susceptibility and 
muon-spin relaxation measurements, finding features at 15 and 65 mT agreeing with the predicted 
boundaries of the magnetization plateau. This suggests that Tb2Ti2 07 is well described by a 
quantum spin ice model with an effective exchange constant of Jeg = 0.17(1) K. 

PACS numbers: 76.75.+i, 75.40.Gb, 75.50.Ee 



Geometrically frustrated magnets are materials in 
which geometrically-induced competition between inter- 
actions prevents local magnetic moments from ordering 
down to temperatures well below the energy scale of those 
interactions. If dynamic magnetic fluctuations persist to 
the lowest experimentally accessible temperatures, the 
low temperature magnetic state is considered a cooper- 
ative paramagnet or spin liquid. These systems provide 
considerable opportunities for experimental tests of the- 
oretical approaches to exotic collective phenomena 
In this context the pyrochlore magnet Tb2Ti207 has 
proved an intriguing conundrum since no magnetic or- 
dering is observed down to the lowest measured temper- 
ature, ^ 50 mK, far below the Curie- Weiss temperature, 
©cw = —19 K 0, and on this basis it has been iden- 
tified as a potential three-dimensional spin liquid [ll-0|- 
Early theoretical models describing Tb2Ti2 07 as an an- 
tiferromagnetic (111) Ising systeiripredicted long-range 
magnetic ordering at around 1 K [J @] , in contradiction 
with the experimental results. One route to overcoming 
this theoretical impasse is to consider the effects of the 
quantum fluctuations of the Tb'^+ magnetic moments, 
which renormalize the low-energy effective Hamiltonian 
of the system from an unfrustrated (111) Ising antiferro- 
magnet to a frustrated (111) Ising ferromagnet the 
quantum analog of the spin ice systems Dy2Ti207 and 
Ho2Ti207 0, Eli- For Tb2Ti207, this quantum spin 
ice model predicts a partial magnetization plateau for a 
magnetic field < 0.1 T applied along the [111] crystal 
axis evident at low-temperatures ^ 0.1 K [12l ]. 

Perpendicular to the [111] axis of a pyrochlore sys- 
tem lie alternating kagome and triangular layers. Apply- 
ing a magnetic field along this axis in the spin ice com- 
pounds Dy2Ti207 and Ho2Ti207 induces a kagome ice 
state with an associated partial magnetization plateau 
that has been investigated in detafl [ii-[l3]. This is anal- 
ogous to the situation predicted for Tb2Ti207, although 



the energy scales are significantly different due to the 
effect of the quantum fluctuations. 

Perturbations are known to drive Tb2Ti207 into a 
magnetically ordered state. Neutron scattering mea- 
surements using flelds applied alon g th e [111] and [110] 
axes 17H21 and applied pressure [22'] found magnetic 
Bragg peaks above 0.125 T and 2 GPa, respectively. The 
magnetic fields used in these studies were, however, too 
large to test the subsequent prediction of a partial magne- 
tization plateau. Measurements applying both pressure 
and magnetic field show that the induced magnetic order- 
ing can be tuned using both parameters [2^. Magnetoe- 
lastic effects also play their part in the low-temperature 
magnetic state of Tb2Ti207. A giant magnetostriction 
]AZ/^] ~ 10~^ is observed at 4.2 K and can be understood 
in terms of the crystal field levels of the Tb^^^ ions [13] . 
X-ray diffraction studies point to structural fluctuations 
in zero magnetic fleld [255 ^^'^ fi^^d ^ cubic-tetragonal 
structural phase transition that can be resolved in mag- 
netic fields S > 25 T [i^. Inelastic neutron scattering 
measurements made in magnetic fields applied along the 
[110] direction show a feature at 0.04 THz indicating the 
presence of a tetragonal lattice distortion [27[- 

To test for the [111] magnetization plateau in 
Tb2Ti207 we have carried out ac susceptibility and muon 
spin relaxation (^SR) experiments using magnetic fields 
applied along the Tb2Ti207 [111] crystal axis. These 
provide complementary information on the magnetic sus- 
ceptibility at two distinct timescales. 

Our ac susceptibility measurements were performed us- 
ing a custom made coil set, thermally anchored to the 
mixing chamber of an Oxford Instruments dilution re- 
frigerator through immersion in liquid ^He. For our /iSR 
experiments [28] the field was parallel to the initial muon 
spin polarization [longitudinal field (LP)] in the ISIS mea- 
surements (0.04 < T < 10 K) and with the initial muon 
spin polarization partially rotated [LP and transverse 
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FIG. 1: (Color online) ac susceptibility data with the field ap- 
plied II [111]: (a) Real part x', (b) Imaginary part x" ■ (Inset) 
Theoretical magnetization curves predicted by the quantum 
spin ice model of Ref. [H for J = 0.167 K. 



field (TF)] for measurements (0.025 < T < 0.9 K) using 
the Low Temperature Fridge (LTF) spectrometer (Paul 
Scherrer Institute, Switzerland). The measured parame- 
ter is the time-dependent muon decay asymmetry, A{t), 
recorded in positron detectors on opposite sides of the 
sample. Having subtracted the background and normal- 
ized the signal, this provides a measure of the spin polar- 
ization P{t) = [A{t) - Abg]/[^(0) - Abg] of the muon en- 
semble as a function of time. For all our fiSK experiments 
the Tb2Ti207 crystals, grown using a floating zone fur- 
nace, were arranged in a mosaic and attached to a silver 
backing plate using a thin layer of GE varnish. The sil- 
ver plate gives a temperature and field-independent back- 
ground signal that can easily be identified and subtracted 
from the asymmetry data. 

The ac susceptibility data are shown in Figure [T] (a) for 
the real part x' ^-nd (b) for the imaginary part x" ■ The 
data recorded at 100 and 125 mK are equivalent within 
error (125 mK data not shown) but there is a clear sepa- 
ration between these data and those recorded at 68 mK 
in both components of the susceptibility. In x' the sepa- 
ration is clear below ~ 200 mT and grows towards zero 
field. This can be com par ed to the theoretically predicted 
magnetization curves |12| in the inset to Fig. [U where a 
separation between the predictions for 50 and 100 mK 
emerges below ~ 300 mT. The temperature dependence 
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FIG. 2: (Color online) (a) Longitudinal field fiSK data at 
50 mK. (b) Transverse field /xSR data at 25 mK. The phase 
shift is due to the orientation of the positron detectors with 
respect to the initial muon polarization. The solid lines are 
the fits to the data described in the text. 



is more pronounced in x" , where at 68 mK the suscepti- 
bility rises to a peak at 15 mT but at > 100 mK it falls 
monotonically. As with x', at 68 mK x" remains distinct 
from that at 100 mK up to ~ 300 mT, in accordance with 
the behavior predicted for the magnetization. Any fea- 
ture at the upper boundary of the magnetization plateau 
is indistinct, also as predicted. Previous measurements 
of the ac susceptibility indicated a partial spin freezing at 
low temperature and a frequency dependence consistent 
with our results Q,!!!]. 

The ^SR data recorded in longitudinal field are shown 
in Fig. [3] (a) and can all be described by the equation: 



Pz(t)=e 



-At 



(1) 



where A is the muon spin relaxation rate. The back- 
ground is weakly relaxing in zero applied field but ef- 
fectively constant in longitudinal applied field, consis- 
tent with the expected behaviour of the silver sample 
holder. That this equation describes the data over the 
whole measured field range demonstrates that the mag- 
netic fields remain dynamic on the timescale probed by 
muons. The transverse field data shown in Fig.[5](b) take 
a similar form except for the muon spin precession. They 
can be described by the equation: 



A^{t) = AsC sin(wsi) -I- Abge 



-Ah 



i(^bgO, (2) 



where the first term describes the signal from the sample 
and the second term describes the background signal. 
The oscillations are sinusoidal rather than the conven- 
tional cosinusoidal behavior because of the detector ge- 
ometry used and the angular precession frequencies are 
related to the magnetic fields experienced by the muons 
as w = (7^/27r = 135.5 MEzT"!). 

The magnetic fluctuations in Tb2Ti207 have previ- 
ously been shown to be fast compared with the range of 
timescales probed by muons [sl. [30l - l32l | . that is to say that 
the system is paramagnetic on the muon timescale, and 
the relaxation rate A can therefore be related to the dis- 
tribution width of magnetic fields at the muon stopping 
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FIG. 3: (Color online) Relaxation rate A measured as a func- 
tion of longitudinal magnetic field applied along the [111] axis. 
Data recorded at 0.5 and 0.725 K are similar to those shown 
for T < 0.2 K so are omitted for clarity. The solid line repre- 
sents the trend for a polycrystalline sample at 0.1 K reported 
in Ref. [sll. The dashed and dotted lines represent fits de- 
scribed in the text. 



site A, the fluctuation time r, and the applied longitudi- 
nal field Blf by the sum of Redfield's equation 33| and 
a field-independent relaxation rate, Aq: 



A 



+ Aq. 



(3) 



We include the field- independent relaxation rate, Aq, be- 
cause our data appear to tend to a field- independent 
value above 0.25 T, which is consistent with the be- 
haviour seen in that field range in Ref. (30| . 

For the T < 2 K longitudinal field data, three regions 
can be identified in the field dependent relaxation rates 
shown in Fig. [31 At small fields up to Bi ^ 15 mT, 
A increases to a peak, then falls steeply to a kink at 
B2 ^ 60 mT, followed by a more gradual fall between 60 
and 250 mT. The fields at which the peak and kink are 
observed are consistent with the boundaries of the mag- 
netization plateau predicted in Ref. [13] (see the inset to 
Fig.[T]). For comparison, we plot the trend A(i?LF) oc 
previously found for polycrystalline data at 100 mK [sij 
as a solid line in Fig. |3l This is similar to the behaviour 
seen in our data above ~ 70 mT. 

Using equation [31 we can test predictions for the field 
dependence of A and r. The simplest assumption is that 
neither depends on field, which effectively describes both 
the 5 and 10 K data. Below 5 K this model does not 
work over the whole field range. However, it is effective 
between the peak {Bi) and kink {B2) shown in Fig. [31 
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FIG. 4: (Color online) Field dependence of the local field 
distribution A estimated from the A values in Fig. [3] using 
equation [l] r = 20 ns, and the Aq values determined from fits 
to the high-field data. 



thereby suggesting that a plateau exists in the local mag- 
netic field distribution in the anticipated field range. Fit- 
ting (with Ao = 0) leads to A ~ 12 mT and r '--^ 13 ns 
at low-temperature, the fit for 25 mK being shown as 
the dashed line in Figure [31 Including the Aq values esti- 
mated using the data above 150 mT (~ 0.5 MHz) reduces 
the value of A by around 15 % and increases t by around 
10 %. Independent of this, the quality of the fits is poorer 
for the 1.5 and 2 K data, which is consistent with ther- 
mal fluctuations breaking up the low temperature state. 
Above ~ 65 mT it is not possible to describe the data 
using equation [31 and the same parameters as between 15 
and 60 mT. This strongly suggests that either the field 
distribution or the fluctuation timescale has changed. 

To visualize the field-dependence of the local magnetic 
field distribution A under the assumption that the fluc- 
tuation timescale r is independent of the weak magnetic 
fields being applied along [111], we can rearrange equa- 
tion [31 into the form: 



y^(A-Ao)(l+7f,SV^)/27^- 



(4) 



Fig.Hlshows the resulting A values assuming r = 20 ns at 
all temperatures. This r value is suggested by the fits to 
the higher temperature data although the form of these 
results is only weakly dependent on r for 12 < r < 25 ns, 
with the plateau region consistently evident in the T < 
2 K data. We expect the fits to equation [31 in the plateau 
region will underestimate r since A is assumed constant, 
which is unlikely at non-zero temperature. The Aq val- 
ues for T < 2 K were estimated by fitting equation [31 
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with constant A and r to the highest field data avail- 
able, fioH ^ 150 niT, and seeking a value that did not 
change significantly with the low-field boundary of the 
fitting window. In the inset to Fig. [1] we show the theo- 
retical magnetization curves from Ref. [l^. There is good 
agreement in the position of the plateau region. Since 
we measure the local magnetic field distribution width A 
rather than M, it is not possible to make quantitative 
comparisons of the magnitude. 

We also carried out muon-spin rotation measurements 
from 2.5 to 250 mT at 25 mK with the aim of finding 
evidence for the level crossing predicted to occur at the 
upper boundary of the magnetization plateau by the the- 
ory described in Ref. [13] ■ No evidence for such an effect 
was found, which could be due to correlations between 
neighbouring tetrahedra or, alternatively, that the pri- 
mary relaxation mechanism for the muon spin is always 
the large distribution of fluctuating local fields. Instead 
we observed a significant negative frequency shift be- 
tween the precession frequencies in the sample and the 
background, K = (wg — ujhg)/^hg- Above ^ 7.5 mT, 
—0.65 < K < —0.625 with negligible field dependence. 
This value is consistent with previous measurements at 
higher fields Sill, which attributed the effect to the 
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large sample magnetization. 

In conclusion, both ac susceptibility and /iSR measure- 
ments of the dynamic magnetic fields in Tb2Ti207 show 
features suggesting that a partial magnetization plateau 
emerges for fields between 12.5 and 65 mT applied along 
the [111] axis. Comparing these fields to the predictions 
of the quantum spin ice model suggests an effective ex- 
change constant of Jeff = 0.17(1) K. The two probes 
show the features being smeared away by thermal fiuc- 
tuations at distinct temperatures - for ac susceptibility 
this is around 100 mK, similar to that predicted for the 
bulk magnetization [12|, whereas for /iSR the features 
persist to around 2 K, consistent with the temperature 
dependence of the relaxation rate observed at a constant 
field in previous measurements |3|. The higher temper- 
ature scale associated with the /iSR results implies that 
on shorter timescales, estimated from fitting the field de- 
pendence in the plateau region to be r ~ 20 ns, the 
fluctuations persist to higher temperature. Such fluctua- 
tions would be quasistatic on the neutron timescale and 
therefore extending previous neutron scattering studies 
to lower flelds would be worthwhile, together with bulk 
measurements of the magnetization. Very recently an al- 
ternative theoretical approach to describing the observed 
properties of Tb2Ti207 was proposed, based on a Jahn- 
Teller-like distortion and a two-singlet system coupled by 
exchange [3J] . Our results provide a further quantitative 
test for this model. 

Parts of this work were performed at the ISIS Facility, 
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tute, Villigen, CH. We thank Andrew Boothroyd for help 
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